Spin and charge ordering in self-doped Mott insulators 
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We have investigated possible spin and charge ordered states in 3d transition-metal oxides with 
small or negative charge-transfer energy, which can be regarded as self-doped Mott insulators, using 
Hartree-Fock calculations on d-p-type lattice models. It was found that an antiferromagnetic state 
with charge ordering in oxygen 2p orbitals is favored for relatively large charge-transfer energy and 
may be relevant for PrNiOs and NdNiOa. On the other hand, an antiferromagnetic state with charge 
ordering in transition-metal 3d orbitals tends to be stable for highly negative charge-transfer energy 
and can be stabilized by the breathing- type lattice distortion; this is probably realized in YNiOa. 



The electronic structure of 'id transition metal oxides 
is described by Zaanen-Sawatzky-AUen (ZSA) scheme j|] 
in which they are classified into two regimes according 
to the relative magnitude of the oxygen-to-metal charge- 
transfer energy A and the d — d Coulomb interaction 
energy U. While the magnitude of the band gap is given 
by U in the Mott-Hubbard regime, it is given by A in the 
charge-transfer regime A < J7. Zd transition-metal ox- 
ides with high valence generally have very small or nega- 
tive charge-transfer energy A and fall in a region which is 
not included in the ZSA scheme Actually, perovskite- 
type 3d transition- metal oxides such as LaCu'^+Oa, ||] 
PrNi^+Og D and SrFe^+Os |] have been studied by 
high energy spectroscopy and have been found to have 
very small or even negative charge-transfer energy A. 

With small or negative A, the highest part of the oxy- 
gen 2p bands can overlap with the lowest part of the up- 
per Hubbard band constructed from the transition-metal 

orbitals so that some holes are transferred from the 3d 
orbitals to the 2p orbitals in the ground state. This state 
can be viewed as a self-doped state of a Mott insulator 
1^ such as has recently been suggested for Cr02 . The 
properties of such system are far from clear and can be 
very rich. It can be a paramagnetic metal, a ferromag- 
netic (FM) metal, and a non-magnetic insulator similar 
to Kondo- insulators [||. However, there exists another 
possibility which has not been explored until now: it may 
have charge ordering or charge density wave. It is possi- 
ble that, in a self-doped state of a Mott insulator, holes 
in the oxygen 2p orbitals undergo charge ordering just 
like doped Mott insulators such as La2-2;Srj;Ni04 In 
this letter, we study this possibility using model Hartree- 
Fock (HF) calculation and show that spin and charge or- 
dered states may appear in perovskites with negative A. 
Based on the calculations, we argue that this phenomena 
occurs in perovskites containing Fe^+ (CaFeOa) and 
Ni3+ (iJNiOa where i? is a rare earth) 0. Specifically, 
we consider the latter system, properties of which, es- 
pecially its strange magnetic properties remain a puzzle 
until now 



We use the multi-band d-p model with 16 Ni sites in 
which full degeneracy of the Ni id orbitals and the oxygen 
2p orbitals are taken into account . The Hamiltonian 
is given by 

H = Hp + Hd + Hpd, (1) 

k,l,(T A;,/>r,cr 

Hd = f^d^d+^^di^ma- + U^dj'^^^d^^rn^d^^^^di^ml 
i,'m(T 2,m 

+ E Kra]di,m^dl^,^d,^„,>y 

+ - /) E d+„^A.rnadl^,J^,m'a 

Hpd= E ykUdt^.P^M+H.c. (4) 

dfrna ^^"6 Creation operators for the 3d electrons at site 
i. d^^^ and p^i„ are creation operators for Bloch elec- 
trons with wave vector k which are constructed from the 
m-th component of the 3d orbitals and from the Z-th 
component of the 2p orbitals, respectively. The intra- 
atomic Coulomb interaction between the 3d electrons is 
expressed using Kanamori parameters, u, u' , j and j' 
]l3[ . The transfer integrals between the transition- metal 
3d and oxygen 2p orbitals V/ff„ are given in terms of 
Slater-Koster parameters (pda) and (pdir). The transfer 
integrals between the oxygen 2p orbitals V^^^' 
pressed by (ppcr) and (ppn). Here, the ratio (pda) / (pdn) 
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is -2.16. (ppcr) and {ppir) are fixed at -0.60 and 0.15, re- 
spectively, for the undistorted lattice. When the lattice 
is distorted, the transfer integrals are scaled using Harri- 
son's law. The charge-transfer energy A is defined by 
e^ — £p + nil, where and ep are the energies of the bare 
3d and 2p orbitals and U {= u — 20/9j) is the multiplet- 
averaged d — d Coulomb interaction. A, U, and (pda) for 
PrNiOa are 1.0, 7.0, and -1.8 eV, respectively, which are 
taken from the photoemission study 

The formally (low -spin d^) compounds i?Ni03 

exhibit a metal-insulator transition as a function of tem- 
perature and the size of the R ion Among them, 
PrNiOs and NdNiOa are antiferromagnetic insulators be- 
low the metal-insulator transition temperature. Neu- 
tron diffraction study of PrNiOs and NdNiOs has shown 
that the magnetic structure has a propagation vector of 
(1/2,0,1/2) with respect to the orthorhombic unit cell 
or is a up- up- down- down stacking of the ferromagnetic 
planes along the (l,l,l)-direction of the pseudocubic lat- 
tice (see Fig. |l|(a)) ||ll[]. In order to explain the magnetic 
structure, orbital ordering of the — and 3z'^ — 
orbitals has been proposed because one of the Cn orbitals 
is occupied in the low-spin configuration pT| ]. How- 
ever, previous model HF calculations have shown that 
the orbital ordered state of — /?>z^ — r^-type has a 
relatively high energy, suggesting that orbital ordering is 
not responsible for the magnetic structure [Q. 

The photoemission study of PrNiOs has shown that 
the charge-transfer energy A of PrNiOa is ~ 1 eV and 
that the ground state is a mixture of the c?^ and d^L con- 
figurations, where L denotes a hole at the oxygen 2p or- 
bitals. Since the ground state has a large amount of oxy- 
gen 2p holes, it is also possible to describe it starting from 
the d^L state. In this picture, the system can be viewed 
as a self-doped Mott insulator and the antiferromagnetic 
and insulating state in PrNiOa and NdNiOa may be in- 
terpreted as a spin and charge ordered state in the self- 
doped Mott insulator. Indeed, our calculations confirmed 
the existence of such ordered states which are consistent 
with the neutron diffraction measurement. They are il- 
lustrated in Fig. 1^. In the state shown in Fig. |l|(a), half 
of the oxygen sites have more holes than the other half. 
The excess holes located at the oxygen sites cause the 
ferromagnetic coupling between the neighboring two Ni 
spins. Therefore, the up-up-down-down stacking of the 
ferromagnetic planes along the (l,l,l)-direction is real- 
ized without orbital ordering. On the other hand, all the 
Ni sites have the same number of "id electrons. Let us de- 
note this state as an oxygen-site charge-ordered (OCO) 
state. In the state shown in Fig. |l|(b), while all of the 
oxygen sites are occupied by the same amount of holes, 
half of the Ni sites have more 3d electrons than the other 
half. This state can be called a metal-site charge-ordered 
(MCO) state. 

In Fig. |^(a), the energies of the spin and charge or- 
dered states relative to the FM and metallic state are 



plotted as functions of the charge-transfer energy A for 
the cubic perovskite lattice. For A < 1 eV, the OCO and 
MCO states exist as stable solutions. The OCO state is 
lower in energy than the MCO state for -5 eV < A < 
1 eV. At A = - 7eV, the OCO and MCO states are al- 
most degenerate in energy. This result indicates that, as 
the charge-transfer energy A decreases, the MCO state 
becomes favored compared to the OCO state. In Fig. 
||(b), the energies of the OCO and MCO states relative 
to the FM state are plotted for the perovskite lattice with 
the orthorhombic distortion which is due to the tilting of 
the NiOg octahedra. Here, the tilting angle is 15° which 
is a typical value found in i?Ni03 ||l^. At A = - 7eV, 
the MCO state is slightly lower in energy than the OCO 
state, indicating that the orthorhombic distortion or the 
GdFeOs-type distortion favors the MCO state. However, 
for A > -5 eV, the OCO state has lower energy than the 
MCO state even with the substantial distortion. Since, 
in PrNiOs and NdNiOa, every Ni site has the same mag- 
nitude of the magnetic moment ||ll[], it is reasonable to 
attribute the antiferromagnetic and insulating state in 
PrNiOs and NdNiOs to the OCO state. In the present 
model calculation without lattice distortions, the OCO 
state is higher in energy than the FM and metallic state 
for realistic A. However, since the charge ordering at 
the oxygen sites is expected to strongly couple with a 
lattice relaxation, a structural modulation may stabilize 
the OCO state as discussed in the following paragraphs. 

The number of 3<i holes N^^ and spin Sd at the Ni sites 
are plotted as functions of A in Fig. ^(a). In the OCO 
state. Nil uniform at all the Ni sites. As A decreases, 
N^J[ becomes smaller because the transfer of holes from 
the oxygen sites to the Ni sites increases. In these solu- 
tions, are approximately two and the population of 
the dx^-y^ orbital is the same as that of the d3z2_^2 or- 
bital, indicating that, in the OCO state, Ni is essentially 
+2 and the orbital degeneracy is lifted. On the other 
hand, are 2.00 and 2.28 in the MCO state for A = 1 
eV. The Ni sites with of 2.00 have the spin Sd of 0.80 
and are Ni2+ as those in the OCO state. The Ni sites 
with N!} of 2.28 have no spin and is well described by 
the d^L_ state which can hybridize with the low-spin d^ 
state. In this sense, the Ni sites can be viewed as Ni'*+- 
like (low-spin d^) sites. Therefore, the MCO state is a 
kind of charge disproportionated state in which two Ni'^+ 
sites are turned into the Ni^+-like and Ni^+-like sites as 
pointed out by Solovyev et al. based on LDA-I-J7 calcu- 
lation Jist . An antiferromagnetic ordering of magnetic 
Ni^+-like sites (see Fig. |l](b)) is also consistent with the 
neutron diffraction results ||ll|. As A decreases, the dif- 
ference of N'J between the Ni'^^-like and Ni'*"'"-like sites 
becomes smaller in the MCO state. The difference al- 
most disappears at A = -7 eV, where the MCO state is 
almost degenerate in energy with the OCO state. Here, it 
is interesting to note that the charge disproportionation 
of 2Fe''+ Fe^+ -|- Fe^+ has been observed in CaFeOa 
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Ipl which has highly negative charge-transfer energy A 
|^|. The number of 2p holes Np at the oxygen sites are 
plotted as functions of A in Fig. ||(b). The OCO state 
has the hole-rich and hole-poor oxygen sites. For A of 1 
eV, Np is ~ 0.33 at the hole-rich oxygen sites and is ~ 
0.22 at the hole-poor oxygen sites. On the other hand, 
in the MCO state, Np is uniform at all the oxygen sites. 

Recently, Medarde et al. have observed strong ^^O- 
"'^^O isotope effect on the metal-insulator transition of 
i?Ni03, indicating that the electron-lattice coupling 
is important in i?Ni03. Very recently, Alonso et al. per- 
formed neutron diffraction studies of YNiOs and found 
the breathing-type distortion which may be an indica- 
tion of charge ordering (it) . In Fig. ^ the relative ener- 
gies of the OCO and MCO states compared to the fer- 
romagnetic state are plotted for A of -1 eV as functions 
of the various lattice distortions with which the charge 
orderings are expected to couple. The MCO state be- 
comes stable with the breathing-type lattice distortion 
as shown in Fig. ^(a). Here, 6o is the shift of the oxy- 
gen ions which gives the breathing-type distortion. The 
MCO state becomes the lowest in energy for rather small 
distortion, indicating that the MCO state coupled with 
the breathing- type distortion is relevant for YNiOa . The 
OCO state can be stabilized with the modulation of the 
bond length which is a consequence of the shift of the 
Ni ions as shown in Fig. |4|(b). Here, the shifts of the 
Ni ion are along the (1,1,1) direction and are given by 
((JNiji^Niji^Ni) and (-(^Nir'^Nir'^Ni)- Consequently, the Ni-0 
bond length for the FM coupling becomes shorter and 
that for the AFM coupling becomes longer. Fig. ^(c) 
shows that the OCO state is also stabilized by the mod- 
ulation of the bond angle which is derived from the tilt- 
ing of the octahedra and the shift of the oxygen ions. 
In this model distortion, the Ni-O-Ni bond angle for the 
FM coupling is 180 ° and that for the AFM coupling is 
smaller than 180 ° . In Fig. §(c), the relative energy is 
plotted as a function of the smaller Ni-O-Ni band an- 
gle. Although these distortions can stabilize the OCO 
state, we need unreasonably large modulations in order 
to make the OCO state lower in energy than the FM 
state. We need more experimental and theoretical in- 
vestigations to identify the lattice distortion realized in 
PrNiOs and NdNiOg. 

In conclusion, we have studied spin and charge ordered 
states in self-doped Mott insulators with small or nega- 
tive charge-transfer energy. It was found that two types 
of charge ordered states are possible: the OCO state with 
charge ordering at the oxygen sites and the MCO state 
with charge ordering at the transition-metal sites. The 
present HF calculation without distortion has shown that 
the OCO state has lower energy than the MCO state for 
moderately small A and that the OCO and MCO states 
are almost degenerate for highly negative A. Since, in 
PrNiOs and NdNiOs, every Ni site has the same mag- 
nitude of the magnetic moment the antiferromag- 



netic and insulating state in PrNiOa and NdNiOs can 
be attributed to the OCO state of the self-doped Mott 
insulators. The OCO state in the self-doped Mott in- 
sulators is novel in that, even without explicit doping, 
the spin ordering at the transition-metal sites and the 
charge ordering at the oxygen sites coexist and couple 
with each other just like the spin and charge ordered 
states in the doped Mott insulators. On the other hand, 
for YNiOa, the strong breathing-type distortion stabi- 
lizes the MCO state. Here, it is interesting to note that 
A of CaFeOa is highly negative and can have the MCO 
state even without the strong beathing-type distortion. 
The charge disproportionated state observed in CaFeOs 
may be regarded as a kind of MCO state in the self- 
doped Mott insulators. In i?Ni03 and CaFeOa, the ho- 
mogeneous state corresponds to an orbitally degenerate 
state {tlgcl for and t^^ej for Fe''+). The charge dis- 
proportionation observed in CaFeOa and the OCO and 
MCO states for iJNiOa may be another way to get rid 
of this orbital degeneracy besides the usual cooperative 
Jahn- Teller (or orbital) ordering. 
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(a) oxygen-site charge ordering (OCO) 



(b) melal-site cliarge ordering (MCO) 



FIG. 1. Schematic spin and charge orderings (a) for 
the oxygen-site charge-ordered (OCO) and (b) for metal-site 
charge-ordered (MCO) states. The arrows and the circles in- 
dicate the Ni spins and the oxygen 2p holes, respectively. The 
larger and smaller circles are for the hole-rich and hole-poor 
oxygen sites, respectively. 



(a) RNiOg without tilting 
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(b) RNiOg witti tilting 
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FIG. 2. Energies of the OCO and MCO states relative 
to the ferromagnetic and metallic state as functions of the 
charge-transfer energy A (a) for the cubic perovskite lattice 
and (b) for the perovskite-lattice with the orthorhombic dis- 
tortion. 



(a) RNiOg witliout tilting 
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(b) RNiOj without tilting 
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FIG. 3. (a) Number of 3d holes A^^ and spin Sd at the 
Ni sites as functions of the charge-transfer energy A. (b) 
Number of 2p holes at the oxygen sites as functions of the 
charge-transfer energy A. 
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FIG. 4. Energies of the OCO and MCO states relative to 
the FM state as functions of (a) the breathing-type distortion 
given by the shift of the oxygen ions, (b) the modulation of 
the bond length driven by the shift of the Ni ions, and (c) 
the modulation of the bond angle driven by the shift of the 
oxygen ions. 5o and ^Ni are the shifts of the oxygen and 
Ni ions. dNi-o is the Ni-O bond length for the undistorted 
lattice. 
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